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Abstract: The related issues of Majorana CP violation in the lepton sector and lepto- 
genesis are investigated in detail in two rather generic supersymmetric models with type I 
see-saw mechanism of neutrino mass generation and A4 flavour symmetry, which naturally 
lead at leading order to tri-bimaximal neutrino mixing. The neutrino sector in this class 
of models is described at leading order by just two real parameters and one phase. This 
leads, in particular, to significant low energy constraints on the Majorana phases 021 and 
a^i in the PMNS matrix, which play the role of leptogenesis CP violating parameters in 
the generation of the baryon asymmetry of the Universe. We find that it is possible to 
generate the correct size and sign of the baryon asymmetry in both A4 models. The sign 
of the baryon asymmetry is directly related to the signs of sin 021 and/or sin 031. 
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1. Introduction 

The presence of two large and one small mixing angles in the lepton sector [1,2], 

sin2 0i2 = 0.304+°:°^^ , sin^ ^23 = 0.501°;^^ , sin^ ^13 < 0.056 (Sex) , (1.1) 

suggests a pattern of neutrino mixing which is remarkably similar to the so called "tri- 
bimaximal" (TB) one [3]. Indeed, in the case of TB mixing, the solar and atmospheric 
neutrino mixing angles 9i2 and ^23 have values very close to, or coinciding with, the best 
fit ones in eq. ( |l.lD , determined in global analyses of neutrino oscillation data, 

sin2 0i2 = l/3, sin2 023 = l/2 . (1.2) 

The TB mixing scheme predicts also that ^13 = 0. Correspondingly, the neutrino mixing 
matrix is given by 

= Utb diag (1, e^"2i/2^ e*"3i/2^ 3) 

where 

/ ,7273 1/^/3 \ 
Utb= -1/^/6 1/^3 I/V2 (1.4) 

\-i/V6 i/Vs -1/V2J 
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and 021 and 031 are Majorana CP violating phases [4,5]. 

TB mixing can appear naturally in models using the tetrahedral group A4 as flavor 
symmetry [6]. ^ The three generations of left-handed leptons and right-handed (RH) 
neutrinos are unified into a triplet representation of the A4 group, whereas the right-handed 
charged leptons are A4-singlets. To be concrete, in the following we focus on two models 
which represent a class of A4 models, namely the models of Altarelli-Feruglio (AF) [8] and 
Altarelli-Meloni (AM) [9]. Both models are based on the Standard Model (SM) gauge 
symmetry group and are supersymmetric. Additional degrees of freedom, flavons, are 
introduced in order to appropriately break the A4 flavor symmetry at high energies. Both 
models have in common that they predict at leading order (LO) a diagonal mass matrix for 
charged leptons and lead to exact TB mixing in the neutrino sector. The mass matrix of 
the RH neutrinos contains only two complex parameters X, Z. Light neutrino masses are 
generated through the type I see-saw mechanism. All low energy observables arc expressed 
through only three independent quantities: a = \3Z/X\, the relative phase 4> between X 
and Z, and the absolute scale of the light neutrino masses. The latter is a combination of 
the unique neutrino Yukawa coupling and \X\ which determines the scale of RH neutrino 
masses. The main diff'erence at LO between the AF and AM models is in the generation of 
the charged lepton mass hierarchy: in [8] an additional Froggatt-Nielsen (FN) symmetry 
U{1)fn is invoked, whereas in [9] the hierarchy between the masses of the charged leptons 
arises through multi-flavon insertions ^. As a result the next-to-leading (NLO) corrections 
in these models are different. The expansion parameter in the A4 models of interest is the 
vacuum expectation value (VEV) of a generic flavon field divided by the cut-off scale A of 
the theory. Its typical size is ~ 0.04 with Ac ~ 0.22 being of the size of the Cabibbo 
angle. 

The fact that the properties of light as well as heavy neutrinos are essentially flxcd by 
the three parameters a, (j) and \X\ leads to strong constraints. We will be interested, in 
particular, in the dependence of the Majorana CP violating phases in the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) matrix [10], which are relevant in leptogenesis and neutrinoless 
double beta decay, on a and (p. Actually, a and 4>, as was shown in [9, 11], are related 
through the ratio r = AmQ/|Am^| = (m^ ~ m,f)/\m,'^ ~ m^j, where Am,Q and Am^ are 
neutrino mass squared differences driving the solar and the dominant atmospheric neutrino 
oscillations. As a consequence, the Majorana CP violating phases depend effectively only 
on a. The related constraints on the neutrino mass spectrum have been studied in [9,11,12]. 

In the present article we investigate in detail the generation of baryon asymmetry of 
the Universe Yb within the AF and AM models. Although our analysis is done for these 
two specific models, it has generic features which are common to models based on A4 
flavour symmetry. In the class of models under discussion the neutrino masses arise via the 
type I see-saw mechanism. Correspondingly, one can implement the leptogenesis scenario 
of matter-antimatter asymmetry generation. As is well known, in this scenario the baryon 

^It can also be derived from models with an 54 flavor symmetry [7]. 

^This is due to a different choice in the AM model of the vacuum alignment of the flavon, (pr, slightly 
different transformation properties of the right-handed charged leptons under A4 and the presence of a 
different additional cyclic symmetry. 
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asymmetry is produced through the out of equihbrium CP violating decays of the RH 
neutrinos (and their SUSY partners - RH sneutrinos P?) in the early Universe [13,14]. 
The observed value of Yb to be reproduced, given as the ratio between the net baryon 
number density and entropy density, reads [15] 

= (8.77 ± 0.24) X 10-^^ (1.5) 



where the subscript "0" refers to the present epoch. 

As has been discussed in [16], the CP asymmetries ej, originating in the decays of 
the RH neutrinos and sneutrinos and and relevant for the generation of the baryon 
asymmetry of the Universe, vanish at LO in the class of models under discussion. Thus, 
Cj 7^ are generated due to the NLO corrections. These are, as already mentioned, different 
in the AF and AM models, so that the results for differ in the two models. 

We calculate the CP asymmetries in the AF and AM types of models and discuss the 
dependence of Yb on the parameter a. This is done in versions of the two models in which 
the RH neutrinos have masses in the range Mj ~ (10^^ 10^^) GeV. As discussed in [9], 
the natural range of RH neutrino masses in the AF and AM models is Mj ~ (10^^ 10^^) 
GeV. Such large values of the RH neutrino masses can lead [17] in SUSY theories with 
see-saw mechanism to conflict with the existing stringent experimental upper limits on the 
rates of lepton flavour violating (LFV) decays and reactions ^, like /x — ^ e + 7, ^ ^ 3e, 
li~ + {A, Z) ^ e~ + {A,Z), etc. The investigation of the LFV processes in the class of 
AF and AM models is beyond the scope of the present study. However, in order to avoid 
potential problems related to the issue of LFV, we work in versions of the AF and AM 
models in which the scale of the RH neutrino masses is lower than in the original AF and 
AM models. This is achieved by minimally extending the AF and AM models through an 
additional Z2 symmetry, which enables us to appropriately suppress the neutrino Yukawa 
couplings. This in turn allows to lower the scale of RH neutrino masses down to (lO^^-hlO^^) 
GeV. 

We perform the analysis of baryon asymmetry generation in the so-called "one flavour" 
approximation. The latter is valid as long as the masses of the RH neutrinos satisfy [18, 19] 
Mj > 5 X 10^-'^(1 -I- tan^ (3) GeV, where tan (3 is the ratio of the vacuum expectation values 
of the two Higgs doublets present in the SUSY extensions of the Standard Model. The 
"one flavour" approximation condition is satisfied for, e.g. tan /? ~ 3 and Afj ~ 10^^ 
GeV. Actually in the models we consider relatively small values of tan/3 are preferable [8]. 
Further, with masses of the RH neutrinos in the range of (10^^ -7- 10^^) GeV one can safely 
neglect the effects of the AL = 2 wash-out processes in leptogenesis [20] . This allows us to 
use simple analytic approximations in the calculation of the relevant efficiency factors rju. 
We perform the calculation of the baryon asymmetry for the two types of light neutrino 
mass spectrum - with normal and inverted ordering. Both types of spectrum are allowed 
in the class of models considered. 

^The indicated problem typically arises if the SUSY particles have masses in the range of few to several 
hundred GeV, accessible to the LHC experiments [17]. 



Yb = 



riB 



n 



B 
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We find that it is possible to generate the correct size and sign of the baryon asym- 
metry Yb in the versions of both the AF and AM models we discuss. The sign of Yb is 
uniquely determined by the sign of sin 0, since all other factors in Yb have a definite sign. 
Interestingly, in the low energy observables only cos ^ is present, so that both sin <^ ^ are 
compatible with the low energy data. Conversely speaking, taking into account the sign 
of the baryon asymmetry Yb we are able to fix uniquely also the sign of sin^, which is 
otherwise undetermined through the low energy data. 

Leptogenesis is not studied for the first time in the class of models of interest. However, 
our work overlaps little with the already existing publications on the subject. In [21] the CP 
asymmetries and Yb were also calculated. This is done, however, not within the context 
of a self-consistent model since instead of computing the NLO corrections, the authors 
introduce ad hoc random perturbations in the Dirac neutrino mass matrix. In [16] only 
the CP asymmetries arc calculated within the AF model without discussing the washout 
effects which can change the results for Yb- In [9] results for the CP asymmetries are 
also given, but the washout effects are not taken into account. Finally, in [22] a highly 
degenerate spectrum of masses of RH neutrinos is considered and the baryon asymmetry 
is produced via resonant leptogenesis. 

The paper is organized as follows: in Section 2 we give a short introduction to the AF 
and the AM models and discuss the changes due to adding the Z2 symmetry. We also give 
the expressions for the NLO corrections relevant for the calculation of the baryon asym- 
metry Yb- In Section 3 we discuss the light and heavy Majorana neutrino mass spectra. 
We study the Majorana phases and their dependence on the parameter a. Our analysis of 
leptogenesis in both models for light neutrino mass spectrum with normal ordering (NO) 
and inverted ordering (10) is given in Section 4. We summarize the results of the present 
work in Section 5. The two appendices contain details on the flavon superpotential and the 
generation of an appropriate VEV for the additional flavon field ^, present in the models 
considered by us. 

2. Variants of the Two Models 

In this section we recapitulate the main features of the AF [8] and the AM model [9]. 
We supplement them with an additional Z2 symmetry to appropriately suppress the Dirac 
Yukawa couplings of the neutrinos and to lower the mass scale of the RH neutrinos. We 
explicitly check that changes in the models connected to the Z2 extension do not affect the 
LO results for the lepton masses and mixings and also only slightly affect the NLO results. 

For an introduction to the group theory of A4 we refer to [8, 9], whose choice of 
generators for the A4 representations we follow. 

2.1 Altarelli-Feruglio Type Model 

In this model the flavor symmetry is accompanied by the cyclic group Z3 and the 
Froggatt-Nielsen symmetry U{1)fn. We add, as mentioned, a further Z2 symmetry to 
suppress the Dirac couplings of the neutrinos. By assuming that the RH neutrinos acquire 
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Table 1: Particle content of the AF variant. Here we display the transformation properties of lepton superfields, 
Minimal Supersymmetric SM (MSSM) Higgs and flavons under the flavor group A4 x Z3 x Z2. I denotes the three 
lepton doublets, e'^, fi'^ and r'^ are the three SU{2)l singlets and v'^ are the three RH neutrinos forming an A4-triplet. 
Apart from v"^ and the flavon Q all fields are neutral under the additional Z2 symmetry. Note that lo is the third 
root of unity, i.e. lo = e~s~ . To accommodate the charged lepton mass hierarchy the existence a U(1)pn, under 
which only the RH charged leptons are charged, is assumed. The U{1)fn is broken by an FN field only charged 
under U{1)fn with charge -1. Additionally, the model contains a C/(l)j{ symmetry relevant for the alignment of 
the vacuum. 



a sign under Z2, the renormalizable coupling becomes forbidden ^. To allow a Yukawa 
coupling for neutrinos at all we introduce a new flavon ^ which only transforms under Z2. 
We call the VEV of ^ = z in the following and assume that z ~ A^A as all other flavon 
VEVs. Clearly, the Majorana mass terms of the RH neutrinos remain untouched, at LO. 
The symmetries and particle content of the AF variant are as given in table 1. At LO 
neutrino masses are generated by the terms^ 



(2.1) 



with (• • • ) denoting the contraction to an ^4-invariant. Thus, the LO terms in the neutrino 
sector are the same as in the original AF model, apart from the suppression of the Dirac 
coupling. Also the LO result that the charged lepton mass matrix is diagonal, is not 
changed compared to the original model. The mass matrices of the neutrinos are of the 
form {niD is given in the right- left basis) 

1 \ ( au + 2bvs —bvs —bvs \ 



■mo = Vu 



1 
1 



-^Vu and ruM 



\ 



—bvs 2bvs au — bvs 
—bvs — bvs 2bvs 



with Vu = {hu), (C) = u and {ifsi) = vs according to the alignment in 
neutrino mass matrix is given by the type I see-saw term 



(2.2) 



The light 



TO. 



(2.3) 



and has the generic size \lvl/K. At the same time, the effective dimension-5 operator 
Ihulhu/A., ^ which can also contribute to the light neutrino masses, is only invariant under 



^Alternatively, one could also let hu instead of transform under the Z2 symmetry to forbid the Dirac 
Yukawa coupling at the renormalizable level. 

^The field ^ does not have a VEV at LO and thus is not relevant at this level. 

^We make the "conservative" assumption that all non-renormalizable operators are suppressed by the 

same cutoff scale A. 
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the flavor group, if it involves two flavons of the type (fg and ^ (^). Thus, its contribution 
to the light neutrinos masses scales as X^v^/K, which is always subdominant compared to 
the type I see-saw contribution. The size of the contribution from the effective dimension-5 
operator is actually of the same size as possible NLO corrections to the type I see-saw 
term. 

Considering the NLO corrections note that these involve for the Dirac neutrino mass 
matrix either the two fiavon combination (^tC the shift of the vacuum of C,. The first 
type gives rise to two different terms 

VAiu^lhs^TKC/^^ + yB{iy^l)3^^ThuC/A-'' (2.4) 

with (• • • )3g(^) standing for the (anti-) symmetric triplet of the product u'^l. The correction 
due to the shift in {() can be simply absorbed into a redefinition of the coupling y,^. Thus, 
using the alignment of tpr as given in [8], the structure of the NLO corrections to is 
the same as in the original model 



SniD = 



f 2yA _ 

-yA - VB ■ (2-5) 

V -VA + VB 



The NLO corrections to the Majorana mass matrix of the RH neutrinos are exactly the 
same as in the original AF model, i.e. 

a S^iu^iy") + d d^ii^^iy") + b (i^^i^^Sips) (2.6) 

+ XA {iy''iy^){ips'PT)/^ + XB {v^V^yiifS^T)" + Xc {v^v")" {ipsVT)' / ^ 

+ XD {iy''u%g{<^s'PT)3s/^ + XE {u''u'')3g{ips<fThA/^ + (l^V'^)3g<^T?/A + Xq {u''u'')3g(pT^/ A 

where Sips, indicate the shifted vacua of the flavons (ps, C and f. Taking 

into account the possibility of absorbing these corrections partly into the LO result, they 

give rise to four independent additional contributions to niM which can be effectively 
parametrized as 

^ 2XD XA XB-xc^ 

SmM = XA XB + 2xc —XD A^A . (2.7) 
\XB- xc -XD XA J 

Compared to these, NLO corrections involving the new fiavon C, are suppressed, since 
invariance under the Z2 symmetry requires always an even number of ^ fields and invariance 
under the Z3 at least one field of the type ips, C or ^. The NLO corrections to the charged 
lepton masses are also the same as in the original model and effects involving ^ can only 
arise at the level of three flavons. 

As we show in appendix A, the VEV of the fiavon ( is naturally also of the order A^A 
as the VEVs of the other flavons and the shift of its VEV is of the size SYEY ~ A^VEV. 
We also calculate its effect on the vacuum alignment of the other flavons and show that 
the results achieved in the original AF model, especially the alignment at LO, remain 
unchanged. 
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Table 2: Particle content of the AM variant. Transformation properties of lepton superfields, MSSM Higgs and 
fiavons under the flavor group A4 X Z4 X Z2 are shown. The nomenclature is as in tabic 1. Apart from hu and the 
flavon C, all fields are neutral under the additional Z2 symmetry. Apart from A4 X Z4 X Z2, the model also contains 
a C/(l)fl symmetry relevant for the alignment of the vacuum, similar to the AF variant. 



2.2 Altarelli-Meloni Type Model 

The AM model, proposed in [9], possesses as flavor symmetry x Z4. To this we add a Z2 
symmetry under which only the Higgs field hu and the new flavon C, transform. Compared 
to the original model, we change the transformation properties of /i„ into 1" under and 
it transforms now trivially under Z4. The flavon is a 1' under ^4 and acquires a phase i 
under Z4. The transformation properties of leptonic superfields, MSSM Higgs and fiavons 
can be found in table 2. The Dirac neutrino coupling is at LO given by 

y,{uH)KC/l^ , (2.8) 

which leads to the same Dirac mass matrix m£, as in the original model, suppressed by the 
factor for z/ K ~ A^. The Majorana mass terms for the RH neutrinos remain unaffected 
by the changes of the model, at LO, 

M{u^u^) + ai{u^v^) + b {vV<fs) (2.9) 

so that the contribution from the type I see-saw to the light neutrino masses arises from 



/lOO\ ( M + au + 2hvs -bvs -bvs \ 



■mo = y, 



V 



Z 



1 I — and txim = \ —bvs 2bvs M + au — bvs 



y 1 y Y -bvs M + au- bvs '^bvs J 

(2.10) 

Thereby the flavon alignment given in [9] is used. Equation ( |2.10D leads to exact TB mixing 
in the neutrino sector. For M ~ A^A, as argued in [9], we find the generic size of the light 
neutrino masses to be A^Vy/A. The effective dimension-5 operator Ihulhu/A arises in our 
variant only at the two flavon level 

i^T^T)"{UK/A' + {^m)'mhl/K^ + {^TVT){U)"hl/A' + {{^T^T)3s{llhs)"hl/A' 
Hi'fml/A^ + i'{^Tll)'hl/K' + e{U)hl/k^ (2.11) 

where we omit order one couplings. Thus, its contributions to the light neutrino masses 
are of order A^v^/A, i.e. of the same size as the expected NLO corrections to the type I 
see-saw contribution, and hence subdominant. 
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The effect of the introduction of the Z2 symmetry and the new field C, on the charged 
lepton sector is the following: an insertion of three flavons, two of which are C,, gives a new 
LO contribution to the electron mass 



(2.12) 



Using the same alignment as in [9], its contribution resembles the one from the operator 
with ^ instead of C, and thus gives also a non- vanishing term in the (11) entry of the 
charged lepton mass matrix. The latter is of the same size as those already encountered 
in the original version of the AM model. Therefore in the variant of the AM model we 
are considering the charged lepton mass matrix is also diagonal at LO and the correct 
hierarchy among the charged lepton masses is predicted. 
At NLO, the Dirac couplings of the neutrinos are 



y,{yH)KKI^ + VAiy'mKlk^ + yB{y''l)zs^si^KI + yc(^^'03^^sC^«/A' ■ (2.13) 



The first two contributions can be absorbed into the LO coupling yi,. Compared to the 
original model, the other corrections are of the same type and generate the same structure 



Note that actually the contribution associated to the coupling is still compatible with 
TB mixing so that only yc can lead to deviations from the TB mixing pattern. As we will 
see in section 4.2, for this reason also the CP asymmetries only depend on the coupling yc- 

All effects to the Majorana mass matrix of the RH neutrinos involving C are negligible, 
since we always need at minimum two fields C and additionally have to balance the Z4 
charge of the operator. Thus, the NLO corrections are only those already present in the 
original model 



The first four contributions can be absorbed into the LO result (or vanish). We do not 
mention effects from shifts in the vacua of cind ^, since these effects can in this model 
also be absorbed into the LO result. The new structures at NLO lead to bmu of the form 



For the charged leptons, additional NLO corrections to the muon and the electron mass 
arise from three and four fiavon insertions, respectively, involving the field C- The operator 



' 2yB -yB - yc -VB + yc \ ^ ^ 
Smr, = -yB + yc "^yB -yB - yc ■ 
V -yB - yc -yB + yc '^yB ) 



(2.14) 



(2.15) 




(2.16) 



3 



(2.17) 
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corrects the muon mass. This type of subleading contribution aheady exists in the original 
model such that no new structures are introduced. The NLO corrections to the electron 
mass are induced through the operator ('^{e'^lipTYhd/A^^, if the shifts in the vacua are taken 
into account, as well as through the four flavon operators 

(2.18) 

All structures arising from these corrections are already generated by the NLO corrections 
present in the original model so that the analysis given in [9] for the NLO corrections is 
still valid in the constructed variant. 

In appendix B we discuss how to give a VEV of the desired size to the field (, the shift 
of this VEV from NLO corrections as well as the effects of C on the flavon superpotential 
of the original model at LO and NLO. 



3. Neutrino Masses and CP Violating Phases in the A4 Models 



The models discussed in the previous section have in common that the Majorana mass 
matrix of RH neutrinos is of the form 



rriM 



(X + 2Z -Z -Z \ 
-Z 2Z X-Z 
-Z X-Z 2Z 



(3.1) 



and the neutrino Dirac mass matrix reads 



/l 0^ 



niD 



1 
1 



A 



(3.2) 



The symmetry of this class of models implies that, at leading order, the neutrino part of 
the Lagrangian depends only on few parameters: X,Z and yy. These parameters are, in 
general, complex numbers. One can set real by performing a global phase transformation 
of the lepton doublet fields. As we will see, CP violating phases, which enter in the CP 
asymmetries of the RH neutrino decays, are functions of the relative phase between X and 
Z. The see-saw mechanism for the neutrino mass generation implies that the full parameter 
space of the neutrino sector can be constrained significantly by the low energy data. 



The RH neutrino mass matrix (3T) is diagonalized by an orthogonal matrix Utb, 
given in ( [1.41) : 

diag(Mle*'^^M2e^'^^M3e*'^3) = U^b^mUtb , (3.3) 



where 



Ml = \X + 3Z| 
M2 = \X\, ip2 



= \X\ |1 + ae*'^| 
arg(A) 



\X - 3Z\ 



\X\ 11 



ae 



i(p I 



"Pi 



arg(X + 3Z) 
:arg(3Z- A). 



(3.4) 
(3.5) 
(3.6) 
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Here a = \3Z/X\ and cj) = arg(Z) - arg(X). 

A light neutrino Majorana mass term is generated after electroweak symmetry breaking 
via the type I see-saw mechanism: 



where 



= — mj^rrijj niD = U* diag (mi, 1712,1^3)11^ (3-7) 



U = iUTB diag (e^'^l/^ e^-^^/^^ e^-^^/^) (3.8) 
and mi^2,3 are the light neutrino masses, 

m.^H^(T)\ -1,2,3 (3.9) 



Mi VA 

The i in eq. ( |3.8D correspond to an unphysical common phase and we will ignore it in what 
follows. We observe also that one of the phases (pk, say ipi, can be considered as a common 
phase of the neutrino mixing matrix, and thus has no physical relevance. In the following 
we always set tpi = 0. 

The parameters a and (f) defined in ( |3.4D -(^^), which determine the RH neutrino 
mass matrix (|3.l| ), can be constrained by the neutrino oscillation data. More specifically, 
we have for the ratio [9]: 

Am| _ (1 + — 2 a cos (?!)) (a + 2 cos </)) /o in\ 

^ ~ \AmW ~ 4| cos(/)| ' 

where AuIq = Am^i = m2 — > and |Am^| = lAmg-^l = |Am|2l are the i/— mass 
squared differences responsible respectively for solar and atmospheric neutrino oscillations. 
Since the value of r is fixed by the data, this relation implies a strong correlation between 
the values of the parameters a and cos4>. Let us note that the sign of sin(/> cannot be 
constrained by the low energy data. As we will see later, the sign of sin cp is fixed by the 
sign of the baryon asymmetry of the Universe, computed in the leptogenesis scenario. 
At 3(7, the following experimental constraints must be satisfied [1]: 

Am| > 

\Am\\ = (2.41 ±0.34) x 10"^ eV^ (3.11) 
r = 0.032 ± 0.006 . 

In Fig. I we show the correlation between a and cos (j), following from ( |3.10| ) taking 
into account ( 3.11| ). Depending on the sign of coscj), the parameter space is divided into 



two physically distinctive parts: cos (p > corresponds to light neutrino mass spectrum 
with normal ordering (NO), whereas for cos^ < one obtains neutrino mass spectrum 
with inverted ordering (10) . 

The main difference between the models we are discussing and the original models 
reported in [8] and [9] is in the mass scale of the RH neutrino fields. In the models 
considered here the predicted RH neutrino masses are always rescaled by the additional 
factor (A^)^ ~ 10~^. Depending on the value of the neutrino Yukawa coupling Ui,, in our 
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Figure 1: The correlation between the real parameter a and the phase 0, which appear in the RH 
neutrino Majorana mass matrix. The figure is obtained by using the 3a range of the parameter r 
given in eq. (3.10). See text for details. 



model the lightest RH Majorana neutrino mass can be in the range from (10^^ -r- 10^^) GeV, 
and up to 10^^ GeV for a neutrino Yukawa coupling y,^ ~ 0(1)- 

For neutrino mass spectrum with NO, the RH neutrino masses show approximately 
the following partial hierarchy [16]: Mi ~ 2M2 ~ IOM3. The lightest neutrino mass mi, 
compatible with the experimental constraints given in ( ^.llD , takes values in the interval 
3.8 X 10~^eV < nil ^ 6.9 x 10~^eV. This implies that the light neutrino mass spectrum is 
with partial hierarchy ^. For the sum of the neutrino masses we have: 

6.25 X 10"^ eV < mi + ms + mg < 6.76 x 10"^ eV . (3.12) 

In the case of 10 spectrum, the overall range of variability of the lightest neutrino mass, 
m3, is the following: 0.02 eV < < 0.50 eV, where only the lower bound follows from 
the low energy constraints. The upper bound was chosen by us to be compatible with the 
"conservative" cosmological upper limit on the sum of the neutrino masses [24,25]. Thus, 
the light neutrino mass spectrum can be with partial hierarchy or quasidegenerate. If the 
spectrum is with partial hierarchy (i.e. 0.02 eV < < 0.10 eV), for the RH Majorana 
neutrino masses, to a good approximation, we have: Mi = M2 = M3/3. Quasidegenerate 
light neutrino mass spectrum implies that, up to corrections ~ 0{r), one has Mi = M2 — 
M3. The sum of the light neutrino masses in the case of 10 spectrum is predicted to satisfy: 

mi + m2 + m3 > 0.125 eV. (3.13) 

We give below the expressions for the lightest neutrino mass in the NO and 10 spectrum 
as functions of a and r. Recall that for fixed value of r, all the parameter space and the 
associated phenomenology is characterized by the parameter a. In the numerical examples 

''This was noticed also in [23]. 
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reported in the following, we always use the best fit value of the ratio r: r = 0.032. By 
expanding with respect to r, we get for the square of the lightest neutrino mass: 



mi 



ml 



^""A^ ( 1 + 2^ ^ (1 + 2a2)3 



2,fJ_ _(l+_a2)r_\ 



NO spectrum 



10 spectrum . 



(3.14) 
(3.15) 



For a = 1 the expression for ml reduces to the one obtained in [9]. 

In the class of models we are considering, the three light neutrino masses obey the 
general sum rule (valid for both types of spectrum) [9, 11]: 



1 

mi 



m2 



(3.16) 



This equation implies a strong correlation between the neutrino masses and the Majo- 
rana phases arising from the RH neutrino mass matrix. The Majorana phases are respon- 
sible for CP violation in leptogenesis and therefore we will discuss them in detail in the 
following subsection. 

3.1 The Majorana CP Violating Phases 

In the following, we use the standard parametrization of the PMNS matrix (see, e.g. [26, 
27]): 



C^PMNS 



/ C12C13 S12C13 Si3e^*'^\ 

-S12C23 - Cl2S23'Sl3e*'' C12C23 - Sl2S23'Sl3e"' S23C13 
\ S12S23 - Cl2C23Sl3e*'^ -C12S23 - Sl2C23Sl3e"' C23C13 / 



. 0^21 ■ '-*^31 

diag(l,e*~,e*~) 



(3.17) 



where Cij = cos 9ij, Sij = sin 9ij, 9ij G [0, 7r/2], 6 G [0, 27r] is the Dirac CP violating phase 
and a2i and 031 are the two Majorana CP violating phases, 021, 31 £ [0, 27r]. From the 



see-saw mass formula (3.7) one can read directly the form of the neutrino mixing matrix 
that arises at leading order in perturbation theory. Taking into account the standard 



parametrization ( 3.17 ) and (1^), the PMNS matrix is indeed: 



C/PMNS = diag(l, 1, -1)Utb diag(l, e'^^/^ e'^'^l'') 
From eqs. ( |3T^ ) and (|1|) we identify the "low energy" Majorana phases as 



"31 



V'2 



(3.18) 



(3.19) 
(3.20) 



We remark that, at this order of perturbation theory, the CHOOZ mixing angle, ^13, is 
always zero as a consequence of the TB form of the neutrino mixing matrix, imposed by 
the broken discrete symmetry. 
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Figure 2: The Majorana phases a2i and a^i in the case of a hght neutrino mass spectrum with 
normal ordering. The parameter r is set to its best fit vahie, r = 0.032. The solutions of equations 
(3.21) and (3.22|) shown in the figure correspond to sin^ < 0. See text for details. 



In the models under discussion the Majorana phases 021 and 031 can also be con- 
strained by using the neutrino oscillation data ^ . After some algebraic manipulation, we 
arrive at the following relations between the CP violating phases 021 and 031 and the 
parameters a and (j) of the model: 

a sin (f) 

tanQ2i = --— 7 (3.21) 

1 + Q cos (j) 

a sin 6 

tana3i = 2^ — ^, 3.22 
a"' — 1 



where a and cos0 satisfy eq. ( 3.10| ). 



In the case of NO spectrum, we have (/> = it e, 27r it e, with e < 0.2 and 0.8 ^ a < 1.2 
(see Fig. 1). If e = 0, the two CP violating phases become unphysical. No CP violation is 
possible in this case. As regards the 10 light neutrino mass spectrum, it is easy to show 
that [9] 2cos(/) ~ —a. The correction, 5o(a), which appears in the right-hand side of this 
equation, is given by 

2ar f 2(l + a>\ 

For light neutrino mass spectrum with inverted ordering, the parameter a varies in the 
interval 0.07 ^ a ^ 2, where the lower limit of a comes from the indicative upper bound 
on the absolute neutrino mass scale used by us, m.1^2,3 ^ 0.5 eV. 

In Figs 2 and 3 we show the behavior of the Majorana phases Q21 and 031 as functions 
of Q, for the NO and 10 mass spectrum, respectively. We choose sincj) < in ( ^.21 ) and 



( 3.22 ). As we will see, this is dictated by reproducing correctly the sign of the baryon 



asymmetry. On the other hand, the relative sign of sin 021 and sin 031 is fixed by the 



*Let us recall that the probabilities of oscillations involving the flavour neutrinos do not depend on the 
Majorana phases of the PMNS matrix [4,28]. Thus, the Majorana phases cannot be directly constrained by 
the neutrino oscillation data. 
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Figure 3: The same as in Fig. 2, but for a liglit neutrino mass spectrum witli inverted ordering. 



requirement that the sum rule in eq. ( |3.16| ) is satisfied. In the case of NO spectrum, the 
phase is close to zero. The maximum value of 021 is obtained for a = 1. At a = 1 we 
have approximately 021 = \/r/3 = 0.1. The other Majorana phase 031 can assume large 
CP violating values. The largest | sina3i| is reached for a = 1: at this value of a we have 
sinasi = —1. 

If the light neutrino mass spectrum is with 10, both phases can have large CP violating 
values. We get sin 021 = 1 and sin 031 = —1 for a « \/2 and a = 1, respectively. 

The Majorana phase 021 can be probed, in principle, in the next generation of exper- 
iments searching for neutrinoless double beta decay [29]. Below we give the expression of 
the effective Majorana mass mge predicted in the class of models we are considering in the 
cases of neutrino mass spectrum with normal and inverted ordering [26]: 



rup 



mi cos 9i2 + \J m{ + Aitlq sin ^126 



NO: 



m^ + \Aml \ |cos2 0i2 + e*"2i 



sm^Oul , 10 . 



(3.24) 
(3.25) 



We recall that in the class of models under discussion, sin^ ^12 = 1/3, cos^ ^12 = 2/3 and 
a non-zero value of ^13 arises only due to the NLO corrections. As a consequence, the 
predicted value of 6*13 is relatively small, ^13 ~ ©(A^ ~ 0.04). Thus, the terms ~ sin^ 6*13 
in ruee give a negligible contribution. Further, since the Majorana phase 021 — (see Fig. 
2), the two terms in the expression for rUee in the case of NO spectrum add up. As a 
consequence, we have: 



mi + -,/ml + Aml 



NO. 



(3.26) 



where 3.8 x lO'^ eV < mi < 6.9 x IQ-^ eV. 

We show in Fig. |^ the effective Majorana mass m-ee and the lightest neutrino mass mi 
versus a, for the NO spectrum. In this case, mee takes values in the interval: 6.5x 10~^ eV < 
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Figure 4: The effective Majorana mass niee (blue continuous line) and lightest neutrino mass mi 
(red dashed line) in the case of a light neutrino mass spectrum with normal ordering. In both cases 
Am^ and r are fixed to their best fit values. 

niee ^ 7.5 X 10~^ eV. Similar conclusion has been reached in [23]. In what concerns the 
10 spectrum, the full range of variability of the effective Majorana mass, compatible with 
neutrino oscillation data, is predicted to be 



For 7713 > 0.02 eV, this implies mee ^ 0.018 eV (see also [9]). 
4. Leptogenesis 

In this section we compute the baryon asymmetry within the AF and AM type models 
defined in Section 2. As we have already noticed earlier, leptogenesis cannot be realized if 
we take into account only the leading order contribution to the neutrino superpotential. In 
order to generate a sufficiently large CP asymmetry, higher order corrections to the Dirac 
mass matrix of neutrinos must be taken into account. The RH neutrino mass spectrum 
in this class of models is not strongly hierarchical. Consequently, the standard thermal 
leptogenesis scenario in which the relevant lepton CP violating asymmetry is generated in 
the decays of the lightest RH (s)neutrino only is not applicable and one has to take into 
account the contribution from the out of equilibrium decays of the heavier RH (s)neutrinos. 
The lepton asymmetry thus produced in the decays of all heavy RH (s)neutrinos {Df 
), i = 1,2,3, is converted into a baryon number by sphaleron interactions. The neutrino 
and sneutrino CP asymmetry e^, which are equal for lepton and slepton final states, is the 
following [30]: 




with mg > 0.02 eV , 10 . 



(3.27) 




(4.1) 
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where 



rriD 



rriD 



(4.2) 



is the neutrino Dirac mass matrix in the mass eigenstate basis of RH neutrinos, and rrii, 
i = 1,2,3 are the hght neutrino masses. The matrix U and the masses mi coincide with 
those given in eqs. ( ^ ) and (^), respectively. The loop function f{mi/mj) is defined as 



fix) 



1 



+ log ( 1 + ^ 



(4.3) 



This function depends strongly on the hierarchy of light neutrino masses. It can lead to 
a strong enhancement of the CP asymmetries if the light neutrino masses rrii and mj are 
nearly degenerate. As we have seen earlier, the neutrinos can be quasidegenerate in mass 
in the case of 10 spectrum. In this case we have to a good approximation f{mi/mj) = 
-f{mj/mi). 

We recall that in the case of 10 spectrum, the lightest two heavy Majorana (s)neutrinos, 
^12 (^12 )' have very close masses. However, the conditions for resonant leptogenesis [31] 
are not satisfied in the models under consideration. Indeed, in all the region of the relevant 
parameter space we have 0.2 ^ a ^ 2. Correspondingly, the relative mass difference of the 
two heavy Majorana (s)neutrinos in question is 



Mo - 



Ml 



mi 
m2 



(2 ^ 14) X lO"'' > max 



{m£)rh) 



DA2 



^2 



10 



-5 



(4.4) 



Under the above condition, the CP asymmetries for each (s)neutrino decay can be computed 
in perturbation theory as the interference between the tree level and one loop diagrams 
(see, e.g. [32,33]). 

The general expression for the baryon asymmetry [34], in which each RH (s) neutrino 
gives a non-negligible contribution, can be cast in the following form: 



Yb = 



riB - riB 



-1.48 X 10~ 



(4.5) 



where r]ij is an efficiency factor that accounts for the effects of washout due to the AL = 1 
interactions of i^j and P| of the asymmetry Y/ . , generated in the decays — > khu, h and 
i>f lihu, k hu- They take into account also the decoherence effects on li caused by the 
and (j 7^ i) interactions. We refer in the following discussion only to the lepton number 
densities. The same considerations apply for the interactions involving slepton states. 

The computation of the efficiency factors in the models under discussion is considerably 
simplified [35,36] (see also [37]). This is due to the fact that to leading order, the heavy 
Majorana neutrinos v^, ^3) ^.s can be shown, couple to orthogonal leptonic states. 

As a consequence, the Boltzmann evolutions of the three lepton CP violating asymmetries, 
associated with the indicated three orthogonal leptonic states, are practically independent. 
Taking into account the above considerations, one can compute the total baryon asymmetry 
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as an incoherent sum of the contributions arising from decays of each of the three heavy 
RH neutrinos: 

3 

YB-Y.'^B^, (4.6) 

i=l 

where 

Yb^ = -1.48 X 10-3ei?7ii . (4.7) 

In the class of models considered the RH neutrino mass scale is set below 10^^ GeV, 
preventing possible washout effects from AL = 2 scattering processes. In this case, the 
efficiency factors rja can be expressed only in terms of the washout mass parameters fhi [20]: 

/3.3xlO-3eV / fhi V'^'^Y^ , 

•"'=(^ ff^ H o.55xlO-^ev J J ■ f"'^' 

where 

fh. . . (4.9) 

Here m/) is the neutrino Dirac mass matrix in the basis in which the Majorana mass matrix 
of RH neutrinos is diagonal with real eigenvalues (see eq. (|4.2| ) . 

4.1 Leptogenesis in the Variant of AF Model 

In this Section we compute the baryon asymmetry for the AF type model, in the one-flavor 
leptogenesis regime. 

In the basis in which the RH Majorana neutrino mass term given in (2.2) is diagonal, 
the relevant matrix that enters into the expression of the leptogenesis CP asymmetries 
(|4.1|) is given by 



- 1 (^\^ 2 2 



+ 



2Re(yA) 2^/2e^^ Re(yA) -^e'^Re{yB) \ 

~ fVT 



2V2e-'^ ReivA) -2 Jfe^^V^ Re{y 



yj=e-'^Re{yB) -2^ le'^"^ Re{yB) -2Re(yA) J 



\AJ\A 



(4.10) 



VuVu 



where and yB are the higher order (complex) Yukawa couplings defined in (2.4) and 
(2.5). We can take all flavon VEVs real without loss of generality. 



ei = -^(y) (^) {6f{mi/m2)sma2iRe{yAf + f{mi/m3)sina3iRe{yBf) (4.11) 



The CP asymmetries {k = 1,2,3) can be written in the following way: 

1 ( Z\'^ /VT 

€2 = (3/(m2/mi)sina2iRe(yA)^ + /(m.2/m.3)sin(a2i - "si) Re(yB)^) 

(4.12) 
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^3 = [j^J [^J (2/(^3/7712) sin(a3i - 021) + /(ma/mi) sin 031) Re{yB) (4.13) 

where ruk are the LO neutrmo masses and z/A « vt/A « A^- Thus, in the model under 
consideration we have 

|efc| oc « 6 X 10-^ k = 1,2,3. (4.14) 

This is the order of magnitude we expect for the CP asymmetry if we require successful 
leptogenesis. Depending on the loop factor f{mi/mj) (eq. 4.3) and the values of the 



Majorana phases, the CP asymmetry can be enhanced or suppressed. 

The washout mass parameters, associated to each of the three lepton asymmetries, are 
given by: 

fhi = mi{l + 0{Xl)) (4.15) 
m2 = m2{l + 0{Xl)) (4.16) 
ms = ms{l + 0{Xl)) (4.17) 

We see that the washout mass parameters, to a good approximation, coincide with the 
neutrino masses. 

Results for NO Spectrum 

We study the baryon asymmetry in the region of the parameter space corresponding 
to a neutrino mass spectrum with normal ordering: 0.8 ^ a < 1.2. The lightest RH Majo- 
rana neutrino in this scenario is z^g. The Majorana phases, that provide the requisite CP 
violation for a successful leptogenesis, are solutions of equations (3.21) and (p.22|) corre- 



sponding to sine/) < 0. The dependence on a of each of the two CP violating phases is 
shown in Fig. 2. We recall that only the solutions corresponding to sine/) < give the 
correct sign of the total baryon asymmetry. 

We show in Fig. 5, left panel, the dependence of the baryon asymmetry Yb on the 
parameter a. The individual contributions to Yb from the decays of each of the three 
RH Majorana neutrinos are also shown. The term Yb3, originating from the lightest RH 
neutrino decays, is suppressed by largest washout effects, with respect to Ybi and Yb2 (see 
(I4T7D )- 



The contribution to the total baryon asymmetry given by Ybi shows an interplay 
between two independent terms proportional to i/a and i/b, respectively. These two terms 
have always the same signs and are of the same order of magnitude. The suppression due 
to the Majorana phase 021 ^0.1 of the term proportional to ua is compensated by the 
enhancement due to the loop factor: we find that (}f{mi/m2)/f{mi/ms) = — (8H-20). The 
same considerations apply to Yb2 ■ Now sin a2i and sin(a2i — «3i ) have the same sign and the 
ratio of the corresponding loop factors is approximately 3f (1712/ mi)/ f {1712/1113) = (20-^30). 

In conclusion, in the case of NO light neutrino mass spectrum, each of the two Majo- 
rana phases 021 and Q31, having values within the ranges allowed by neutrino oscillation 
data (see Fig. 2), can provide the CP violation which is required in order to have suc- 
cessful leptogenesis. Even in the case in which the term proportional to sin 031 in the CP 
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Figure 5: AF type model: baryon asymmetry versus a in the cases of neutrino mass spectrum with 
normal (left panel) and inverted (right panel) ordering. In each plot we show: i) the total baryon 
asymmetry Yb (red continuous curve), ii) Ybi (green dashed curve), iii) Yb2 (orange dotted curve) 
and iv) Yes (blue dot-dashed curve). On the right panel, the lines corresponding to Ybi and Yb2 
overlap. In both cases sin < and Ato^ and r are fixed at their best fit values. The results shown 
in the left (right) panel correspond to i/a = 2.5 and = 3 (j/^i = 0.4 and ys = 2). The horizontal 
dashed lines represent the allowed range of the observed value of Yb, Yb G [8.5,9] x 10^^^. 



asymmetries is strongly suppressed (which corresponds to the case of strong fine-tuning of 
ys ^ 1), successful baryogenesis can be naturally realised for values of the Majorana phase 
OL2\ ~ (0.04 0.10) and a moderately large neutrino Yukawa coupling da ~ (2.5 3.0). 

Results for lO Spectrum 

We now study in detail the region of the parameter space for which the neutrino mass 
spectrum is with inverted ordering and is hierarchical. This scenario is realized for 0.2 < 
a < 2. In the following, we report the behavior of the baryon asymmetry in all the interval 
of variability of a, compatible with an 10 neutrino spectrum and for which the computa- 
tion of the CP asymmetry can be done in perturbation theory. Thus, the results we show 
for 0.07 < OL < 0.2 should be valid provided the renormalisation group (RG) effects [38] 
are sufficiently small in the indicated region. 

In Fig. |5|, right panel, we plot the different contributions to the baryon asymmetry, as 
we have done previously for the normal hierarchical mass spectrum. 

The Majorana CP violating phases which enter into the expressions for the CP asym- 
metries are reported in Fig. 3. The solutions of equations ( |3.21| ) and ( 3.22| ) corresponding 



to sin(/> < must be used also in this case in order to obtain the correct sign of the 
baryon asymmetry. Now f| is the heaviest RH Majorana neutrino and the washout ef- 
fects for the CP asymmetry generated in the decays of this state are less strong since 
they are controlled to LO by the lightest neutrino mass 7713: ms = 7713. We note, how- 
ever, that also in this scenario the contribution of the term Ybz in Yb is always much 
smaller than the contribution of the other two terms Yb\ and Yb2- This is a consequence 
of the strong enhancement in the self energy part of the loop function that enters into 
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Figure 6: The same as in Fig. 5 but for the AM type model. The Yukawa coupling yc is set to 
2/(7 2 in both the figures. 

the expressions for Ygi and Yb2- Indeed, if the spectrum is inverted hierarchical, we have 
f{mi/m2) — —f{'m2/mi) ~ 50/(m3, m-i^2)- For this reason the CP violating phase 031 
gives, in general, a subdominant contribution in the CP asymmetries ei and £2 when the 
Yukawa couplings yA and ys are of the same order of magnitude. This conclusion is valid 
even in the region of the parameter space where a^i ~ 37r/2. 

The analysis of all the parameter space defined by q, compatible with low energy 
neutrino oscillation data, in the AF type model, shows that in both the normal and inverted 
patterns of light neutrino masses, the Majorana phases can provide enough CP violation 
in order to have successful leptogenesis, even in the case in which only one of the phases 
021 and 031 , effectively, contributes in the generation of the CP asymmetry. 

4.2 Leptogenesis in the Variant of AM Model 

In this section we study the generation of the baryon asymmetry of the Universe in the 
variant of the AM model considered by us. We work in the one-flavour leptogenesis ap- 
proximation. The quantity relevant for the calculation of the CP asymmetries in this case 
is: 




1 ^ 




(4.18) 



+ 




where the ys and yc are defined in eqs. (2.13) and (2.14). Again we can choose all flavon 
VEVs to be real without loss of generality. 



The CP asymmetries in this model are given by 




(4.19) 
(4.20) 
(4.21) 



where mi^s are again the LO neutrino masses (see eq. (|3.9| )). The leptogenesis CP violating 
phase now coincides with the Majorana phase asi. Moreover, the CP asymmetries ei^s 7^ 
are controlled by only one parameter, yc-, of the matrix of neutrino Yukawa couplings 
(2.14), the reason being that only this parameter breaks the TB form of the latter. 

As we see, the heavy RH Majorana neutrino z^f "decouples": the CP violating lepton 
asymmetry is produced in the out of equilibrium decays of the heavy Majorana neutrinos 
and 1^3 alone. This constitutes a major difference with the variant of the AF model, 
analyzed in the preceding subsection. After the lepton asymmetries are converted into a 
baryon asymmetry by sphaleron processes, the final matter-antimatter asymmetry of the 
Universe can be estimated as: 



where lei, for z = 1,3, are given in eq. ([l.Tp. The LO washout mass parameters mi^s are 
the same as in the variant of the AF model: 



in the cases of neutrino mass spectrum with normal and inverted ordering. Both types 
of spectrum are allowed in the model considered. The ranges of possible values of the 
Majorana phase 031 which provides the correct sign of the baryon asymmetry are shown 
for the NO and 10 spectra in Figs. 2 and 3 (they are the same as for the AF type model). 

We observe that, as in the variant of the AF model, the suppression of the term Yes 
with respect to Ybi in the case of NO spectrum is due to the relatively larger washout effects 
in the generation of the asymmetry . The maximum of the total baryon asymmetry Yg 
is reached for a ~ 1 where the CP violating Majorana phase 031 = 37r/2. (see Fig. 2, right 
panel). 

In what concerns the 10 spectrum, the two terms Ybi and Ybz enter with the same sign 
in the total baryon asymmetry and are of the same order of magnitude. The enhancement 
of the asymmetry for a < 0.7 is explained by the increase of the loop function /(mi/m^) = 
—/(m^/mi) in the region of quasi-degenerate light neutrino mass spectrum. 

In this class of models, successful leptogenesis can be naturally realized for both types 
of spectrum - NO and 10, for an effective Yukawa coupling yc ^ 1.5. 

5. Summary 

In the present work we studied the related issues of Majorana CP violating phases and 
leptogenesis in variants of two prominent (and rather generic) super symmetric A4 models 



Yb = Ybi + Ybs 



(4.22) 



mi = mi{l + 0{Xl)) 



(4.23) 
(4.24) 



In Fig. I 



we show the dependence of the baryon asymmetry on the parameter a 
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[8,9] which naturally lead at leading order (LO) to tri-bimaximal (TB) mixing in the lepton 
sector. The pattern of neutrino mixing suggested by the existing neutrino oscillation data 
is remarkably similar to the TB one. Both models are supersymmetric and employ the type 
I see-saw mechanism of neutrino mass generation. They predict at LO a diagonal mass 
matrix for charged leptons and lead to exact TB mixing in the neutrino sector. The mass 
matrix of the RH neutrinos contains only two complex parameters X, Z . All low energy 
observables are expressed through only three independent quantities: the real parameter 
a = \3Z/X\, the relative phase (p between X and Z, and the absolute scale of the light 
neutrino masses. The latter is a combination of the neutrino Yukawa coupling and the 
parameter \X\ which determines the scale of RH neutrino masses. 

The main difference between the original models and those considered by us is in the 
scale of RH neutrino masses. In the original models this scale is around (10^^ 10^^) GeV. 
In order to avoid possible potential problems with LFV processes we consider versions of 
both models in which the scale of RH neutrino masses is lower, namely, is in the range 
of (10^^ 10^^) GeV. This is achieved by imposing an additional Z2 symmetry capable of 
suppressing sufficiently the neutrino Yukawa couplings. As a consequence, the mass scale 
of the RH neutrinos is lowered as well. We discussed in detail the flavon supcrpotcntial in 
the modified models of interest. The results obtained at leading order and next to leading 
order (NLO) in the original models arc still valid in the extensions wc consider. 

The two Majorana phases of the PMNS matrix, 0:21 and 031, effectively play the role 
of leptogenesis CP violating parameters in the generation of the baryon asymmetry. In 
the models considered both the phases 0:21 and asi and the ratio r = AmQ/|Am^| are 
functions of the two parameters a and (f). We analyzed in detail the dependence of the 
two "low energy" Majorana phases a2i and 0:31 on a and (p. In contrast to the low energy 
observables, like neutrino masses and the effective Majorana mass in neutrinoless double 
beta decay, ruee, we show that these phases depend both on smcp and coscp, and not only 
on cos (f). We show also that the sign of the baryon asymmetry Yb uniquely determines the 
sign of sin cf), which has to be negative: sin cf) < 0. 

In the case of neutrino mass spectrum with normal ordering (NO), 021 is shown to 
be small, 021 ^ 0.1. In the types of models considered sin^ ^13 is also predicted to be 
small, sin^ 6*13 ~ 10^"^. As a consequence, the contributions of the terms oc sin^^i3 in 
niee are strongly suppressed. The lightest neutrino mass is predicted to lie in the interval 
(3.8 6.9) X 10~^ eV, thus the neutrino mass spectrum is with partial hierarchy. The 
effective Majorana mass rriee has a relatively large value, rUee ~ 7 x 10~^ eV. We note that 
if a2i had a value close to tt, one would have rUee <S 10~^ eV. Depending on a, the phase 
a^i can take large CP violating values. For light neutrino mass spectrum with inverted 
ordering (10), the Majorana CP phases 0121 and 031 vary (for sin^ < 0) between and vr 
and TT and 2tt, respectively. 

Throughout this study we have neglected renormalization group effects on neutrino 
masses and mixings which can be large for a quasi-degenerate (QD) light neutrino mass 
spectrum. A QD spectrum can arise in the models considered if Am^ < (i.e. the 
spectrum is with inverted ordering) and a < 0.2. However, this corresponds only to a 
small portion of the parameter space of the models. 
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As has been already discussed in the literature, in the models of interest leptogenesis 
is not possible at LO: the corresponding CP asymmetries Ci vanish. Thus, the inclusion of 
NLO effects is crucial for the generation of the baryon asymmetry Yb ■ More precisely, the 
NLO effects correcting the neutrino Dirac mass matrix itid give rise to non- vanishing 
and therefore to non- vanishing Yb- Due to this fact the CP asymmetries are naturally of 
the order of 10~^ (independent of the precise value of the loop function). Further, although 
the AF and AM type models lead to the same results at LO, they differ at NLO so that 
the CP asymmetries generated in the two models are different. 

We find that it is possible to generate the correct size and sign of the baryon asymmetry 
Yb in the versions of both the AF and AM models we discuss. The study of leptogenesis 
was performed in the framework of the one flavor approximation and by using analytic 
formulae for the relevant efficiency factors rja. Since the mass spectrum of the RH neutrinos 
is generically not strongly hierarchical, the decays of all three RH (s) neutrinos contribute 
to the generation of the baryon asymmetry. We find that the correct magnitude as well as 
the correct sign of the baryon asymmetry Yb can be easily obtained in the AF and AM 
type models for most values of the parameter a and natural values of the NLO couplings. 
As already mentioned, the sign of Yb uniquely fixes the sign of sin (j). The latter cannot be 
determined by low energy observables since they exhibit only a cos 0-dependence. 

To conclude, the results of our detailed study show that SUSY models with flavour 
symmetry and type I see-saw mechanism of neutrino mass generation, which gives rise to 
tri-bimaximal mixing and Majorana CP violation in the lepton sector, can account also 
successfully for the observed baryon asymmetry of the Universe. 
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A. Flavon Superpotential in the AF Type Model 



In the construction of the flavon superpotential we closely follow [8] and introduce an addi- 
tional U{l)ji symmetry under which driving fields have charge +2, superfields containing 
SM fermions +1 and flavons, h^j and FN ficld(s) are uncharged. To give a VEV of order 
A^A to C we introduce a new driving field (^o which is a singlet under all symmetries of 
the model, apart from carrying a U{1)r charge +2. The terms contributing to the flavon 
superpotential containing at LO read ^ 

tf^ = M^Co + 9aCoC^ + gbCoiv^Tfr) ■ (A.l) 
Analogously to the original model, we demand a vanishing F— term of 

+ gaC^ + gbivll + 2ipT2VT3) = . (A.2) 

At the same time, the field (, does not couple to the other driving fields, ~ (3,1), 
~ (3,0;^) and ~ (1)^^) under {A4,Zs), in the model at LO. Thus, their F— terms 
read as in [8]. We find as solution 

= --{M^ + gbV^) (A.3) 

ga 

and the same results for the VEVs of tpr, ffis, ^ and ^ as in [8]. For the mass parameter 
being of order A^A the VEV z is also of order A^A. 

Concerning the NLO contributions stemming from ^ to the alignment of the flavons 
^T, <^5, C and ^ we find just one term 

|CM^W (A.4) 

which gives an additional contribution 

/ /\/f?\ .,2 

(A.5) 

to the shift Sv^i of (px- Its size is A^A, as expected. Furthermore, the shifts Svt2,3 remain 
unchanged and thus still equal. The shifts in the vacuum of (ps and ^ are also unchanged 
and the VEV of ^ is still a free parameter. 
The NLO terms affecting w'^ read 

^^^=i^2^^I' (A-6) 

i=l 

with 

II = CoiPTPTPT) , li = CoifS^SPs) , I3 = CoCiPSPs) , I4 = CoiifSfs) , 

(A.7) 




® Terms such as Qahuhd are not relevant, since we assume that the flavor symmetry is broken much above 
the electroweak scale. 
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The result for the shift in the VEV of (, z + 6z, in the usual linear approximation, is 
^9939a V J zA 2ggi 

with 5(4 = —^4 and = 8^3 as introduced in [8]. This shift 5z in {() is of order A^A. 
Additionally, we find that - unless some non-trivial relation among the couplings in the 
flavon super potential is fulfilled - the VEVs of all driving fields vanish at the minimum. 

B. Flavon Superpotential in the AM Type Model 

In order to induce a VEV for the flavon ( we add a driving fleld Co which transforms as 1' 
under A4, with —1 under Z4 and is invariant under the Z2 symmetry. Since it is responsible 
for the vacuum alignment, its charge under the ^7(1)^ symmetry is +2. The LO potential 
for Co is of the form 

wi = gaCo^ + 9bU'PW)" + 9cUi'? ■ (B.l) 
Prom the i^-term of Co we can derive 

9a(^ + 9h{vli + 2(^Tim) + 9c{i'f = . (B.2) 

Thus, z takes the value 

~ {9bvl + <7c(n')^) = (f# + 9) {u'f (B.3) 

9a 9a \ ^"■2 / 

SO that z (X u' holds in case of no accidental cancellations, u' is a free parameter in [9] 
which is taken to be of order A^A. 

As one can check, the field C does not have renormalizable interactions with the driving 
fields, ipQ ~ (3, —1), ~ (3, 1) and ^0 ~ (1) 1) under {A4, Z4), of the original model. Thus, 
the results for the vacuum alignment found in [9] still hold. 

At NLO the field C contributes to the fiavon superpotential of the original model 
through 

ICi^IvsY , (B.4) 

while it does not introduce any contribution at this level involving (p^ or ^o- 

The NLO effects on the vacuum alignment of the field ( stem from (order one coeffi- 
cients are omitted) 

^CoC'^ + jCoiV^m^s)" + jm'PT^T)" + jCoC' {^T^SY + jCo^'^'C . (B.5) 

Computing the effect of all NLO terms on the vacuum alignment one finds that still all 
shifts Svsi are equal, i.e. the shifts do not change the structure of the vacuum, that the 
generic size of all shifts - for mass parameters and VEVs of order A^A - is A^A and the free 
parameter u' is still undetermined. 

Eventually, we checked that all driving fields can have a vanishing VEV at the mini- 
mum. 



z^ 
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